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The  work  accomplished  by  the  Principal  Investigator  under  grant  AFOSR 
84-0159  from  its  initiation  on  July  1,  1984  through  the  present  may  be 
classified  into  four  broad  areas,  each  of  which  is  discussed  in  turn. 

A.  Statistical  Inference  for  Repairable  Systems. 

In  Whitaker  and  Samanlego  [16],  we  focus  on  the  development  of 
estimation  procedures  for  the  life  distribution  F  of  a  new  system  based  on 
data  on  system  lifetimes  between  consecutive  repairs.  The  'Imperfect 
Repair'  model  introduced  by  Brown  and  Proschan  [5]  postulates  that,  at 
failure,  the  system  is  repaired  to  a  condition  as  good  as  new  with 
probability  p  and  is  otherwise  repaired  to  the  condition  just  prior  to 
failure.  In  treating  Issues  of  statistical  inference  for  this  model,  we 
begin  by  pointing  out  the  lack  of  identif lability  of  the  pair  (p,F)  as  an 
index  of  the  distribution  of  interfailure  times  Ti,T2,...  We  show  that 
data  pairs  (Ti,Zj^),  i  -  1,2,...  are  necessary  to  render  the  parameter  pair 
(p,F)  identifiable,  where  is  a  Bernoulli  varaiable  which  records  the 
mode  of  repair  (perfect  ox  imperfect)  following  the  i^^  failure.  We 

A 

demonstrate  that  the  nonparametric  maximum  likelihood  estimator  of  F 
exists  only  in  special  cases,  but  that  a  neighborhood  maximum  likelihood 
estimator  always  exists  and  may  be  derived  in  closed  form.  We  demonstrate 

A 

the  strong  uniform  consistency  of  F  under  mild  assumptions,  and  prove  weak 

A 

convergence  of  an  appropriately  scaled  version  of  F  to  a  Gaussian  process. 

X 
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These  results  are  shown  to  apply  to  various  experiaencal  designs 
(including  renewal  testing  and  inverse  saapling)  and  to  extend  to  the  age- 
dependent  inperfect  repair  model  of  Block,  Borges  and  Savits  [2]. 

Llndley  and  Singpurwalla  [8]  have  studied  the  modeling  of  failure 
time  data  from  experiments  involving  changing  environments .  They 
consider,  for  example,  the  situation  in  which  two  components  are  tested 
in  a  controlled  environment  and  subsequently  tested  in  a  (perhaps 
harsher)  operating  environment.  Letting  Exp(X)  represent  the  exponential 
distribution  with  density 

f(t)  -  Ae‘^*^,  t>0, 

they  analyze  the  model  in  which  component  i  is  assumed  to  be  distributed 
according  to  Exp(Ai)  in  laboratory  experiments  and  according  to  Exp(f}Xj^) 
under  operating  conditions,  where,  typically,  ij>l.  In  Currit  and 
Singpurwalla  [6],  this  model  is  studied  further,  and  the  mechanics 
necessary  for  Bayesian  estimation  of  A  (assuming  A]^-A2**A)  are  developed. 
Samaniego  [12]  provides  a  complementary  analysis,  focusing  on  the  maximum 
likelihood  estimators  of  the  parameters  of  the  exponential  distributions 
governing  the  lifetime  of  components  observed  under  both  laboratory  and 
operating  conditions.  The  treatment  is  general,  applying  to  any  number  of 
components  observed  under  dual  conditions.  Some  potential  applications  of 
this  analysis  exist  in  modeling  repairable  systems.  For  example,  if  a 
collection  of  parallel  systems  in  k  non* Identical  components  are  tested 
until  failure,  and  if  these  systems  are  repaired  and  retested,  it  seems 
reasonable  to  assume  that  the  pre-  and  post-repair  environments  satisfy 
the  Llndley-Singpurwalla  assumptions. 
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B.  Inference  Besed  on  Observed  Extreae  Values. 

Samaniego  and  Whitaker  [14]  treat  experiments  In  which  only  record- 
breaking  values  (e.g.,  values  smaller  than  all  previous  ones)  are 

observed.  The  data  available  may  be  represented  as  Xi,Ki,X2,K2 . 

where  X]^,X2,...  are  successive  minima  and  K]^,K2>...  are  the  numbers  of 
trials  needed  to  obtain  new  records.  Hie  problem  of  estimating  Che  mean 
of  an  underlying  exponential  distribution  is  considered,  under  both  fixed 
sample  size  and  inverse  sampling  schemes.  Under  inverse  sampling,  we 
demonstrate  certain  globl  optimality  properties  of  an  estimator  based  on 
the  "total  time  on  tes"  statistic.  Under  random  sampling,  it  is  shown 
than  an  analogous  estimator  is  consistent,  but  can  be  improved  for  any 
fixed  sample  size. 

When  only  a  single  sequence  of  random  records  are  available, 
efficient  estimation  of  the  underlying  distribution  F  is  possible  only  in 
a  parametric  framework.  However,  under  the  assumption  that  Che  process 
of  observing  random  records  csan  be  replicated,  Samanlego  and  Whitaker 

A 

[15],  derive  and  study  the  nonparametric  maximum  likelihood  estimator  F 
of  F.  The  strong  uniform  consistency  of  this  estimator  is  established  (as 
the  number  of  replications  grovs  large),  and  its  asymptotic  distribution 

A 

theory  is  described.  The  performance  of  F  is  compared  to  that  of  two 
reasonable  competing  estimators. 

Boyles  and  Samanlego  [3]  have  studied  nonparametric  estimation 
under  nomination  sampling,  a  sampling  process  in  which  every  observation 
is  the  maximum  of  a  random  sample  from  some  population.  Assuming  Chat 
all  samples  are  taken  from  a  single  Whderlylng  distribution  F,  data  may 
be  viewed  as  consisting  of  pairs  (X^.K^),  where  is  the  size  of  the  i^^ 
sample  and,  given  •  kj^,  X^  is  distributed  according  to  F^^.  Willemain 


[18],  who  discussed  nomination  sampling  In  the  context  of  health  care 
delivery,  proposed  an  estimator  for  the  median  of  F  under  the  assumption 
that  ■  N,  a  fixed  Integer.  In  our  context,  the  assumption  of  a  fixed 
sample  size  N  from  each  population  Is  relaxed;  with  K  taken  as  random, 
the  problem  of  nonparametrlc  estimation  of  the  distribution  function  F 
Is  considered.  The  nonparametrlc  maximum  likelihood  estimator  of  F  Is 
the  derived.  Its  consistency  Is  demonstrated,  and  Its  asymptotic  behavior 
a  stochastalc  proacess  Is  Identified.  Conditions  are  given  under  which 
these  asymptotic  results  hold  for  nonrandom  K.  A  by-product  of  this 
development  is  the  consistency  of  tflllemaln's  estimator  of  the  median. 
Several  applications  are  considered.  For  example,  nomination  sampling 
arises  naturally  in  certain  reliability  experiments;  the  applicability  of 
derived  estimator  in  problems  Involving  designed  redundancy  is  noted. 

A  detailed  analysis  of  data  on  34  yearly  maximum  floods  of  the  Nidd 
River  is  presented,  and  the  estimator  of  the  underlying  flood 
distribution  F  is  displayed.  In  Boyles  and  Samanlego  [4] ,  a  "structural 
approach"  to  the  estimation  of  F  is  considered,  and  the  alternative 
method  Is  shown  to  be  especially  well  suited  for  semlparametrlc 
formulations  of  the  problem,  i.e.,  for  situations  in  which  K  is 
appropriately  modeled  by  a  discrete  parametric  family  of  distributions, 
while  F  Is  unrestricted. 

C.  Inference  for  Nonparametrlc  Reliability  Models. 

Much  of  the  recent  literature  In  statistical  reliability  theory 
focuses  on  nonparametrlc  modeling  and  Inference  for  life  distributions. 
The  nonparametrlc  approach  to  modeling  involves  constructing  classes  of 
life  distributions  which  reflect  a  particular  qualitative  aging  behavior 
of  components  or  systems.  The  appeal  of  these  classes  Is  that  a  few 


5 


natural  constraints  suffice  to  capture  intuitive  notions  of  aging.  On  the 
other  hand,  statistical  Inference  results  have  generally  been  more 
difficult  to  obtain  for  these  models.  Ue  have  been  successful  in  treating 
a  number  of  Inference  problems  of  this  type.  Three  specific  problems  are 
described  below. 

(1)  Let  X  be  a  positive  random  variable.  The  distribution  F  of  X  Is 
said  to  be  "new  better  than  used  in  expectation,"  i.e.,  NBUE,  if  E(X)  2: 
E(X-t|X>t)  for  all  t^O.  Suppose  Xj^,...,Xn.  is  a  random  sample  from  an 
NBUE  distribution  F.  Uhltaker  and  Samanlego  [17]  consider  the  problem  of 
estimating  F  by  a  distribution  which  Is  itself  NBUE.  The  estimator  Gji, 
defined  as  the  NBUE  distribution  supported  on  the  sample  which  minimizes 
the  (sup  norm)  distance  between  the  NBUE  class  and  the  empirical 
distribution  function,  is  studied.  The  strong  \inlform  consistency  of 
is  proven,  and  a  numerical  algorithm  for  obtaining  G^  is  given.  ■ 

(ii)  Suppose  the  survival  function  S  of  a  random  variable  X  is  a  member 
of  the  NBU'Cto)  class.  By  definition  (see  Hollander,  Park,  Proschan 
[7]),  S  must  satisfy  the  following  condition: 

(1)  S(x+to)  S  S(x)S(To)  for  all  x  ^  0. 

Let  [Xi,...,Xn}  be  a  random  sample  from  S.  The  goal  is  to  find  an 
estimator  that  itself  belongs  to  the  NBU-Ctg}  class  and  has  the  same 
asymptotic  optimality  properties  as 


Reneau  and  Samaniego  [9]  study  the  estimator  S  of  S,  given 
implicitly  by  Sn(x)  0<x<to 

“  min  [s„(x),  S(x-to)S(to)]  x  >  t,,  ’ 


A 

It  is  shown  that  S  may  be  written  in  teir-  of  the  empirical  survival 


curve  Sn  as  follows: 
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"  OrSk^fx/tg]  Sjj(x-lcto)}  , 

where  [u]  denotes  the  greastes  Integer  less  than  or  equal  to  u. 

A 

Reneau  and  Sananiego  [9]  establish  the  fact  that  the  estimator  S  is  a 
strong  uniformly  consistent  NBU-Cto)  estimator  of  S  which  converges  to  S 
at  the  best  possible  (pointwise  and  mean  square)  rates.  Confidence 

A 

procedures  based  on  S  are  also  discussed.  ■ 

(iii)  For  any  set  A  C  [0,«),  we  define  a  survival  function  S  to  be  New 
Better  Than  Used  with  respect  to  the  set  A  (NBU-A)  if 

S(x|t)  <  S(x|0)  for  all  x  i  0,  t  e  A, 
where  S(x|t)  -  P(X>x+t|X>t) . 

This  class  clearly  contains  the  NBU  class  of  distributions,  and  it  is 
in  turn  contained  in  the  NBU-Ctl  class  for  all  t  e  A.  Reneau  and 
Sananiego  [10]  study  recursive  estimation  formulas  for  a  number  of  special 
cases.  When  A  is  a  finite  set.  the  resulting  estimators  are  shown  to  be 
strong  uniformly  consistent  estimators  of  the  corresponding  survival 
curves.  The  same  is  true  for  the  set  A  ••  CktQ:  k-1,2,...).  When  the  set 
A  contains  a  limit  point,  the  methods  used  on  the  preceeding  cases  are 
not,  in  general,  successful  in  producing  a  consistent  estimator.  Even  in 
these  cases,  however,  they  do  provide  useful  estimators  in  one  important 
area  of  application,  namely,  when  products  are  subjected  to  bum*in,  with 
those  items  which  fall  before  a  predetermined  time  6  discarded.  The 
survival  function  obtained  by  applying  left- truncation  to  a  survival 
function  S, 


S^(x) 


{ 


1 

S(x)  /  S(«) 


is  itself  in  the  NBU-A  class  whenever  S  is. 


if  0  s  x  <  { 

if  £  ^  X 

Thus,  random  sampling  from  a 


NBU-A  distribution  under  fixed  left-truncation  such  as  bum-  in  yields  a 


survival  function  of  the  desired  type.  When  S  has  the  property  that  S(£) 
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-  1  for  some  f  >  0,  the  recursive  approach  to  estimating  S  It  shown  to 
lead  to  a  strong  uniformly  consistent  estimator  with  optimal  convergence 
rates.  ■ 

D.  Structural  Results  in  Reliability. 

We  have  obtained  a  nximber  of  general  results  which  clarify  the 
behavior  of  coherent  systems  and/or  positive  (failure  time)  random 
variables.  Among  these,  the  following  are  worth  special  mention. 

(I)  Samanlego  [11]  derives  a  representation  for  the  failure  rate  of  an 
arbitrary  coherent  system  when  the  lifetimes  of  its  components  are 
Independently  distributed  according  to  a  common  absolutely  continuous 
distribution  F.  The  system  failure  rate  is  written  explicitly  as  a 
function  of  F  and  Its  failure  rate.  The  representation  Is  used  In  several 
examples,  including  an  example  showing  that  the  well-known  closure  theorem 
for  k-out-of-n  systems  in  i.i.d.  IFR  components  proven  by  Barlow  & 

Proschan  can  be  extended,  but  not  to  all  coherent  systems.  The  class  of 
coherent  systems  for  which  such  closure  obtains  is  characterized.  ■ 

(II)  Samanlego  [13]  discusses  various  methods  and  criteria  for  comparing 
coherent  systems.  A  theoretical  result  Is  derived  for  comparing  systems 
of  a  giver,  order  when  components  are  assumed  to  have  Independent  and 
Identically  distributed  lifetimes.  The  result  gives  sufficient 
conditions  for  the  lifetime  of  one  system  to  be  stochastically  larger 
than  that  of  another  system.  ■ 

(ill)  The  identity  E(X)  -  /^(x)dx,  where  X  Is  a  nonnegative  random 
variable  and  F(x)  -  P(X>x)  is  well  known  in  reliability  theory  and 
survival  analysis.  In  a  1982  paper,  Samanlego  extended  this  identity  to 
produce  an  expression  of  the  k*^^  moment  of  a  nonnegative  random  variable 
in  terms  of  an  Integral  involving  the  survival  function; 
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E(X*‘)  -  ...  F(xi+. .  .+xvt)dxi. .  .dxfc, 

when  the  expectation  exists.  This  identity  can  be  used  in  establishing 
certain  inequalities  for  the  New  Better  Than  Used  (NBU)  class  of 
distributions.  Arnold,  Reneau  and  Samanlego  [1]  generalize  the  identity 
above  or  arbitrary  multivariate  random  variables.  They  then  use  the 
result  to  derive  inequalities  for  multivariate  extensions  of  the  family 
of  NBU  distributions. 
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